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We present the results of an Ultracold neutron (UCN) production experiment in a pulsed neutron 
beam line at the Los Alamos Neutron Scattering Center. The experimental apparatus allows for a 
comprehensive set of measurements of UCN production as a function of target temperature, incident 
neutron energy, target volume, and applied magnetic field. However, the low counting statistics 
of the UCN signal expected can be overwhelmed by the large background associated with the 
scattering of the primary cold neutron flux that is required for UCN production. We have developed 
a background subtraction technique that takes advantage of the very different time-of-flight profiles 
between the UCN and the cold neutrons, in the pulsed beam. Using the unique timing structure, we 
can reliably extract the UCN signal. Solid ortho-D2 is used to calibrate UCN transmission through 
the apparatus, which is designed primarily for studies of UCN production in solid O2. In addition 
to setting the overall detection efficiency in the apparatus, UCN production data using solid D2 
suggest that the UCN upscattering cross-section is smaller than previous estimates, indicating the 
deficiency of the incoherent approximation widely used to estimate inelastic cross-sections in the 
thermal and cold regimes. 
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I. INTRODUCTION 

Ultracold Neutrons (UCN) are free neutrons with max- 
imum energy of 300 neV, equivalent to ~ 1 mK. Their 
kinetic energy is so low that they may be contained and 
accumulated in material bottles, magnetic traps, and the 
Earth's gravitational field, for durations up to hundreds 
of seconds The kinetic energy is comparable to the 
Zeeman splitting in a magnetic field of a few Tesla, and 
neutron state with 100% spin polarization can be pre- 
pared through simple field filtering. Storability and easy 
polarization of UCN make them the tool of choice in 
many experiments to measure fundamental properties of 
the neutron to unprecedented precision. 0, Q The most 
precise measurements of the neutron electric dipole mo- 
ment [H and the neutron beta-decay lifetime [5|, @ all 
employ UCN in a trap. Any innovation that increases 
the storable UCN density and the deliverable UCN flux 
will simplify many technical difficulties in these challeng- 
ing experiments. For example, a larger UCN density 
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would reduce the size of most EDM experiments, making 
the implementation of stringent magnetic field uniformity 
and large electric field a less daunting task. In addition 
to applying UCN for the studies of fundamental physics, 
there are also many applications to condensed matter 
physics and possible enhanced sensitivities derived from 
the long wavelength of UCN0|, if a more intense source 
were to become available. 

UCN are already present in the Maxwell-Boltzmann 
spectrum of thermalized neutrons emerging from fission 
reactors. However, the percentage of the low-energy pop- 
ulation is so low that even the most powerful research re- 
actors cannot easily deliver high enough UCN density to 
be of interest. The use of a cold neutron moderator could 
shift the energy distribution to as low as 30 K, resulting 
in an increased UCN flux. In addition, many tricks have 
been implemented, including gravitational deceleration 
[H and Bragg deflection on mechanical turbines to further 
slow down neutrons Q . To increase the phase-space den- 
sity beyond the limit imposed by the Liouville theorem, 
the most efficient method is to dissipate the neutron en- 
ergy through excitations in condensed matter. Phonons 
in many materials have energies (~ meV) comparable to 
that of moderated cold neutrons. Golub and Pendlebury 
[l(J first proposed a "superthermal" source to use super- 
fluid helium as a UCN converter. In this type of source 
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the neutron gives up all of its kinetic energy by exciting 
collective excitations in the interacting material. The de- 
pletion of the UCN population via absorption of the same 
excitation energy (so-called "upscattering" ) can be sup- 
pressed to an arbitrarily small level simply by lowering 
the temperature of the converter. In the superthermal 
source, the number density of UCN could accumulate for 
durations as long as the /3-decay lifetime of free neutrons. 

To preserve a sizable population of UCN (~ mK) in 
a source environment of a few Kelvins, it is essential 
to avoid re-thcrmalization by delaying upscattering. It 
is also important to control sources that could lead to 
UCN loss, such as nuclear absorption. Any medium 
possessing small nuclear absorption and significant neu- 
tron scattering cross-sections, regardless of the details 
of the excitation, are potential candidates. Successful 
implementation of this idea have been demonstrated in 
superfluid helium [ITl - fl3j and solid deuterium p^4l7| . 
Both of these UCN converters have phonon excitations 
that match well with the incident cold neutrons for ef- 
ficient energy transfer through single inelastic scatter- 
ing. Other materials could potentially make even more 
efficient UCN converters. For one, 16 has a neutron 
absorption cross-section five times smaller than that of 
2 H. In addition to phonons, solid oxygen (S-O2) in its 
low temperature phases has strong magnetic interac- 
tions which could be harnessed for UCN production [l8| . 
These magnetic excitations have been widely studied us- 
ing Raman scattering [IH [2(| as well as neutron scatter- 
ing [2T| - [23| | . It has been shown that the low tempera- 
ture phases possess very different dynamics that often 
entangle translational, librational, and spin excitations 
through the orientation-dependent couplings between the 
diatomic oxygen molecules on the solid lattice sites. 

We present an apparatus used to test S-O2 as a UCN 
converter. A background subtraction technique has been 
developed to enable collection of UCN production data in 
a pulsed neutron facility, where the cold neutron back- 
ground can be significant. To extract the absolute ef- 
ficiency of UCN production from the source, we care- 
fully characterize the efficiencies of UCN transmission 
and UCN detection. We calibrate the efficiency of the 
apparatus using the well-studied solid ortho-deuterium 
(0-D2) as a standard reference. Here, we present a com- 
prehensive test of the physics of UCN production in o- 
D2 using this apparatus. The experimental results of 
UCN production in S-O2 (including dependence on in- 
cident neutron energy, source volume, and the applied 
external magnetic fields) will be presented in forthcom- 
ing papers. 



II. EXPERIMENTAL METHODS 

A. Overview 

The experiment is carried out on a cold neutron beam 
line, where the cold neutron flux and energy spectrum 



are well characterized. The neutron flux is collimatcd 
and monitored continuously. Even though the available 
flux is significantly smaller compared to a previous exper- 
iment using a dedicated mini-spallation target (13 , fl5j , 
the removal of the uncertainties associated with the pro- 
duction and transport of the cold neutron flux is essen- 
tial in understanding the results of UCN production. The 
new apparatus consists of a neutron scattering target cell 
which allows us to grow solid target in-situ. The cell is 
placed in a cold neutron beam line to act as the UCN 
production source. 

UCN, once produced inside the source, are extracted 
and then detected in a neutron detector placed suffi- 
ciently far away from the primary neutron beam. The 
strategic placement of the detector helps to control the 
backgrounds. Elastically scattered neutrons from the in- 
tense cold neutron beam are the dominating source of 
background in this experiment. Determining the UCN 
extraction and transport efficiency is the major chal- 
lenge in determining the absolute UCN production cross- 
section. We therefore employ an 0-D2 converter to bench- 
mark the performance of the apparatus. Production of 



UCN in 0-D2 has been demonstrated at PNPI [24[ and 
the sup erthermal behavior was measured at LANSCE 
[Til Il5l |. Further neutron transport and UCN production 
cross-sections in 0-D2 have been carefully carried out at 
PSI [Hill and by the Mainz/Munich group [H . 

A unique feature of the apparatus (illustrated in Fig.rTJ 
is the integration of a high magnetic field, which allows 
for the study of the magnetic excitations in S-O2. The 
target cell is placed in the center of a superconducting 
(SC) solenoid with the axis aligned along the direction 
of the neutron beam. A cylindrical stainless steel guide 
serves as the UCN transport guide, as well as the insu- 
lating vacuum for the cryogenic target cell. The guide 
extends beyond the warm bore of the SC solenoid and 
connects to the rest of the UCN guide system. 

The cell is larger than the cold neutron mean free path 
(MFP) in 0-D2 so that the majority of incident cold neu- 
trons scatter at least once before escaping the cell. Inside 
the filled cell, a small fraction of these scattering events, 
less than 1 part in 10 6 , results in the production of a 
UCN. Some of the produced UCN escape from the cell, 
enter the guide system, and move toward the detector. 

The UCN guide system was made of stainless steel 
tubes 10.16 cm in diameter. The pharmaceutical-grade 
guides are internally electro-polished to about ~ 0.25 /im 
mean variation of the surface roughness. The total length 
was approximately 2 m with two 90-degree bends to di- 
rect the UCN perpendicularly out of the primary cold 
neutron beam, then vertically downward into the detec- 
tor. The target cell is cylindrical, with the cold neutron 
beam impinging on the front surface and the UCN pri- 
marily extracted from the opposite surface downstream. 

Ostensibly, only UCN and very cold neutrons (VCN), 
the latter with a lower probability, could totally reflect 
from the internal walls of the guide system and propagate 
to the detector. However, we found a significant back- 



3 



ground due to cold neutrons which arc diffracted and/or 
scattered incoherently by the section of UCN guide that 
intersects the beam. Even though the probability is quite 
low (~3x 10~ 8 ), the large flux of cold neutrons passing 
through the apparatus leads to a background comparable 
to the UCN count rate. In sec IIV1 we will show that this 
background can be removed by using its very distinct 
timing structure. With the relatively short guide sys- 
tem, the apparatus has a high UCN throughput, without 
being plagued by the cold neutron background. 

Eventual use of S-O2 as a UCN converter requires good 
temperature control of the cryogenic target cell, essential 
to reproducing the thermodynamic conditions for repeat- 
able crystal growth. The major challenge arises from the 
reduction of the molar volume of S-O2 by 12.5% upon 
cooling through the three solid phases (7, /3, a) at satu- 
rated vapor pressure. A discrete 5% change at the 7-/3 
transition [2g, [2?} at 44 K is the bottleneck for attaining 
large-sized S-O2 cryocrystals in the j3 and a phases. Poor 
crystal quality, which adversely affects UCN extraction, 
needs to be controlled in order to gain understanding of 
the physics of UCN production in S-O2. Density fluc- 
tuations (caused by cracks) over a range comparable to 
and/or larger than the UCN wavelength could result in 
additional scattering and lead to a reduced MFP Q • Sev- 
eral groups have observed this effect in S-D2. In partic- 
ular, Atchison and coworkers [l6| measured a ~20 barn 
increase of the total cross-section in 0-D2 after several 
thermal cycles. There are no data on s-Cv 




FIG. 1: A schematic of the apparatus. Cold neutrons (1) en- 
ter the front Aluminum window (2) and are incident on the 
cryogenic target cell (3). Ultra cold neutrons (UCN) produced 
within the cell diffuse out of the cell and into a highly polished 
stainless steel UCN guide (4), which direct them through two 
90 degree bends into a UCN detector (5). The cell itself re- 
sides a 5.5 T (maximum) superconducting solenoid, housed 
inside a liquid helium cryostat (6). The target cell is cooled 
by a pulse tube refrigerator (7). A cut-away view of a partially 
filed cell (3) is shown in the inset. 



Unlike o-D 2 , the MFP of UCN in s-0 2 is not limited 
by the incoherent scattering length, and is theoretically 
infinite because of the zero nuclear spin of the oxygen 
nucleus 8 0. On the other hand, any additional voids and 
cracks could limit the MFP of UCN and significantly alter 
the extraction efficiency. With this apparatus, we have 
investigated different factors limiting the MFP of UCN 
in the production target by changing the length of the 
target cell along the cold neutron beam axis. 



B. Cold Neutron Beam 

The UCN production apparatus was constructed and 
tested at the Indiana University Cyclotron Facility, and 
then installed on the Flight Path 12 (FPI2) in the Lujan 
Center at LANSCE in August 2008. This neutron facil- 
ity generates pulsed neutrons at 20 Hz from a spallation 
target. FPI2 is coupled to a liquid hydrogen neutron 
moderator with straight neutron guides. Details of the 
neutron source emission time distribution, guide p erfor- 
mance, and overall intensity can be found in [28l - l30f . The 
neutron spectrum peaks at 3.3 mcV (^40 K) [3l[ l32l ]. 
The neutron energy is determined from the timc-of-flight 
(TOF), t, over the flight path length, L, by 



E(t) = ^mL 2 r' 2 

with the corresponding energy resolution (5E) of 
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(i) 



(2) 



Here 5t is the emission time, and 5L is the uncertainty 
in the flight path length. The arrival of the proton pulse 
defines the to, however, details of moderator geometry, 
neutron slowing-down and diffusion within the neutron 
source itself limit the timing resolution to the emission 
time [28l | . The emission time is the time spread for neu- 
trons of a given energy. The fastest neutrons observed by 
TOF are around 100 meV, with a resolution of around 
2% for St ~ 150 fis and L =21.1 m. The relative un- 
certainty in the flight path length is negligible. A frame 
overlap chopper absorbs the long wavelength neutrons 
and restricts the lower limit of the incident neutron en- 
ergy to E>1.2 meV. The most recent measurement [30| 
reported the neutron flux at the sample position to be 
(2.0 ± 0.1) x 10 7 n cm _2 s _1 , integrated over the spec- 
trum from 1.2 to 20 meV, for 100 fj,A proton current on 
the spallation target. This measured intensity includes a 
boost in flux from an m=3 guide. 

Typical neutron spectra collected with our instrument 
on FP12 are shown in Fig. [2J using an empty cell and 
a cell filled with o-D 2 at 4.81±0.02 K. The charge- 
integrated voltage signals from the cold neutron monitors 
are digitized with a 12-bit waveform digitizer with a sam- 
ple rate of 102.4 kHz. The neutron intensity I(t) is nor- 
malized to the average proton beam current of 100 fiA. 
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FIG. 2: An example of the neutron monitor signals I(t) 
recorded with the waveform digitizer with an averaged proton 
current of 100 fj,A on the spallation target. The empty cell 
signal (black) is shown together with the 4 K 0-D2 filled cell 
data (red) for both the incident beam monitor (Mi) and the 
transmitted beam monitor (M2.) The solid gray lines show 
the features which are aligned to determine the monitor posi- 
tion relative to the incident beam monitor, which is fixed at 
21.11 m. Vertical axis is the voltage output from the ampli- 
fier, and its scale is set by gain resistor in the amplifier circuit. 
The monitor signal I(t) is the integrated voltage signal per 
time channel normalized to the duration of data-taking, i.e., 
1000 s. 



A few features in the beam line modifies the expected 
Boltzmann spectrum of neutrons from the moderator. 
First, the signal was zero at low energies because the 
frame-overlap chopper was closed to prevent the slower 
neutrons from appearing in the next frame. The neu- 
tron flux increased as the chopper was opened, leading 
to a rising edge in the spectrum. The two most pro- 
nounced dips in the spectrum were the Al [111] and Al 
[200] Bragg edges from aluminum structural elements and 
vacuum windows in the cold neutron beam path. Addi- 
tional peaks due to Bragg scattering from solid D2 and 
the stainless steel UCN guides in our apparatus appeared 
in the downstream monitor (M2). 

The neutron flux monitor Mi was mounted at the end 
of the cold neutron guide, and was kept at a fixed position 
of 21.11±0.03 m downstream from the moderator. By 
aligning the Bragg edges, we calibrated the position of 
M2 to be 22.75±0.04 m from the moderator, with the 
error bar on the distance determined from the width of 
the Bragg peak. Neutron signals measured by Mi were 
unchanged between runs, indicating good stability. 

The flux of transmitted cold neutrons (collected in 
M2) is used to calculate the total scattering cross-section 
of the target material (0-D2) during calibration runs. 
Several factors related to the non-ideal geometry need 
to be taken into account in calculating the total cross- 
sections. To avoid blocking the fill/vent line below solid- 
ification temperature and thus prevent pressure hazards, 
the cell was only partially filled with liquid (volume frac- 
tion / = 0.601 ± 0.014 of o-D 2 ). The fill volume was dc- 
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FIG. 3: The total cross-section of 0-D2 derived from Eq.© 
is compared with measurement done by (35|. The theoretical 
cross-section (solid line) assuming a hexagonal closed packed 
(HCP) structure is shown for comparison. 



termined from the pressure drop in the storage volume, 
together with the known density of the liquid. With the 
cell uniformly illuminated, the neutron flux measured in 
M2 thus consists of both un-attenuatcd neutrons from the 
upper part of the cell and the attenuated neutrons from 
the lower part of the cell. We determine the cross-section 
using the following algorithm: 



na{E)x = In 



J filli 



(1-/) 



(3) 



where a is the total neutron cross-section per molecule, 
n is the molecular number density, and x is the target 
thickness (3.56 cm). The combined effects of digitizer 
voltage resolution and frame overlap chopper restrict our 
measurement in energy range from 1.2 to 30 meV. 

The total cross-section for 0-D2 is shown in Fig. [31 
together with the data of Seiffert [33| . We obtain the best 
agreement with the previous measurements by adjusting 
the fill fraction to / = 0.57 ± 0.02. This reduced value 
reflects the expected 4 % fractional change of volume 
upon cooling 0-D2 to 4 K 34]. However, this volume 
change could also represent a 1.6 ± 0.7 mm thick layer 
of solid frozen to the walls in the unfilled, upper portion 
of the cell. Given the high vapor pressure of 0-D2 , both 
effects likely contribute to the filling fraction / smaller 
than that estimated by the pressure drop record. Finally, 
we note the cell is about 1.6 MFP thick, and thus multiple 
scattering was likely to occur, leading to an excess in the 
measured cross-section as observed in the range below 
10 meV. 
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Gas Handling System and Cryogenic Target 
Cell 



The materials of interest to this experiment, 0-D2 and 
O2 gas, are both highly volatile, and thus much care was 
put into the design and the construction of the Gas Han- 
dling System (GHS). A schematic of GHS is shown in 
Fig. |U The GHS was constructed out of all stainless 
steel VCR tubing and connectors with all-metal seals. 
Only dry pumps and fomblin pumps were connected to 
the GHS. In the beginning of each run, it was cleaned by 
pumping until the pressures were less than 10~ 4 Pa mea- 
sured on the GHS panel before any gas was introduced. 
The cell was filled from a 200 1 storage tank via a flow 
controller on the GHS. A check valve was installed be- 
tween the cell and the storage tank. It was set to open at 
a 1.3 x 10 4 Pa pressure differential to relieve overpressure 
in the cell back to the storage tank to contain the volatile 
gas. In addition, a drop-out plate on the UCN guide sys- 
tem (which contains the cell) mitigated any catastrophic 
overpressure in the event the cell should burst. 

There was only one fill line (0.3175 cm diameter SS 
tubing) that connected the GHS to the cell. In order to 
prevent blockage, the fill line had to be kept at tempera- 
tures above the solidification temperature of the material 
under study. In addition, the fill line was connected to 
the cryogenic cell using a stainless steel flange. The low 
heat conductivity of the stainless steel flange allowed for 
a large temperature gradient, isolating the cryogenic cell 
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TABLE I: Cell geometries used in the experiment (see inset 
in FigdJ 



Cell Type 


Base Length (cm) 


Cup Length (cm) 


Volume (cc) 


Small 


1.143 


(flat window) 


38.82 


Medium 


1.143 


2.413 


117.4 


Large 


6.198 


2.413 


290.2 



FIG. 4: Schematic of the gas handling system (GHS). 



from the warm fill line. The fill line was wrapped with a 
nichrome wire heater and temperature-controlled with a 
PID loop. 

Inside the cell, the solid was frozen from the liquid 
phase. The triple points are 18.7 K, 16.9 kPa for D2 and 
54.4 K, 152 Pa for O2. The 0-D2 gas was borrowed from 
the dedicated UCN source for the UCNA experiment [36[ . 
The ortho-para ratio of the gaseous D2 was measured via 
Raman spectroscopy at the beginning and the end of the 
experimental run. The initial measurements showed a 
contamination of 2.9 ± 0.02 % para-D 2 and 0.17 ± 0.11 % 
hydrogen dcutcride (HD). The para fraction increased 
0.2 % over the 50-day duration of the experiment, and 
the HD impurity level remained constant within error. 
After the D2 calibration run, the D2 gas was thoroughly 
pumped out before the oxygen gas was introduced. O2 
was supplied via a high pressure gas cylinder at 99.999% 
purity. 

The target cell consists of a type 6061 Al cylindrical 
cup (several lengths available) bolted to a circular base 
flange, and installed horizontally in the magnet bore. 
The cell wall downstream of the beam direction con- 
sists of a thin Al window (0.254 mm) to facilitate UCN 
transmission and subsequent detection. To locate the 
cell within the center of the magnetic solenoid, the top 
flange is attached to one end of a high purity (99.999%) 
Al bar ^30 cm in length. The Al bar is cooled by a pulse 
tube refrigerator with a total cooling power of 1.5 W at 4 
K. The target cell is made vacuum-tight with an indium 
wire seal. Throughout the experiment, three different cell 
lengths were used with the base flange/cup combinations 
shown in table [Jj The inner diameter of the cup side of 
the cell is 6.4 cm, and the base flange side is 6.67 cm. We 
used the medium cell for the 0-D2 calibration runs. 

The cell was surrounded by two layers of cold shields. 
The inner cold shield was connected to the cold head 
(cold stage) of the pulse tube refrigerator, while the outer 
shield was cooled by the warm stage (~ 50 K). The whole 
assembly was placed inside the stainless steel UCN guide, 
which runs through the warm bore of a separate helium 
cryostat that houses the superconducting magnet. The 
inner cold shield was lined with a thin nickel foil (50 /im) 
to increase the collection efficiency of UCN emerging from 
the sides of the cell. Several G10 rings were used to center 
the cold shields and the cell, preventing thermal shorts 
between the cold shields and the inner wall of the UCN 
guide. Stable temperature operation with fluctuations no 
larger than 3 mK was achieved from 4.8 to 80 K with an 
empty cell. When filled, fluctuations of each temperature 
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sensor increased to 15-20 mK. In addition, the tempera- 
ture gradient across the cell when filled with solid o-D 2 
increased to ~ 60 mK. The UCN guide was evacuated 
to less than 10~ 4 Pa to provide insulation required to 
operate the cryogenic target cell. Throughout the exper- 
iment, a residual gas analyzer constantly assayed the gas 
composition in the guide vacuum, to monitor for leakage 
from the cell. 

We used Lakeshore CERNOX 1050-SD Resistive Tem- 
perature Detectors (RTDs) for thermometry because of 
their resistance to high radiation and insensitivity to high 
magnetic fields. Temperature readout and heater power 
was controlled using a Lakeshore 330 Temperature Con- 
troller. We developed a MATLAB based DAQ to imple- 
ment slow controls for the RTDs, heaters, and pressure 
sensors. The same DAQ also controlled a fast digitizer 
to read out the cold neutron monitors, and the multi- 
channel scalar (MCS) used to collect TOF spectrum from 
the UCN detector. This DAQ program thus provided a 
centralized platform for thermometry control and moni- 
toring, data acquisition, visualization, as well as the on- 
line data analysis. 



III. RESULTS 




Pulse Height Channel (Volts) 

FIG. 5: The pulse height spectrum measured by the multi- 
wire proportional counter. Solid dots are the measured data. 
The red curve is the theoretical fit of the spectrum including 
the wall effects from the daughter nuclei of proton and triton. 
The solid blue line is the fitted gamma background. With the 
threshold set to enclose only the full energy peak, the UCN 
detecting efficiency is 85 ± 3 %. 



Ultracold Neutron Detection 



MCS. To collect enough data for statistical analysis, TOF 



An improved version of the multi-wire proportional 
counter [37| was used to detect UCN. Extensive neutron 
shielding using 0.5 m thick borated polyethylene and B4C 
powders surrounding the detector was essential to reduce 
the background rate from 300 Hz (un-shicldcd) to 134 ± 
2 mHz. The counter, filled with 1 kPa of 3 He gas, has a 
small, yet non-zero efficiency of 1.3 % to detect cold neu- 
trons with a 40 K Maxwellian spectrum. A typical UCN 
detector count rate with a 60% filled cell of o-D 2 at 5 K 
with 100 fiA proton current is 612 ± 13 mHz. The beam- 
off ambient neutron backgrounds were only 20 mHz, and 
thus most of the 134 mHz observed is due to the trans- 
mitted cold neutrons which were elastically scattered in 
the UCN guide. 

Inside the FP12 cave, there was a high level of 7 ra- 
diation as a result of neutron captures. The 7s were 
detected by the UCN detector as low energy pulses, but 
their intensity was high enough to prevent clean separa- 
tion between the 7 background from the higher energy 
neutron peak. We set the counter threshold to read the 
full energy neutron peak. The large threshold enhanced 
7 rejection to nearly 100%, at a cost of reducing the effi- 
ciency for neutron detection to 85 ± 3 %. Note that the 
detector front window also attenuates the UCN flux, the 
effect of which will be included as a part of the transmis- 
sion efficiency. 

Signals from the UCN detector were recorded using 
a multi-channel scalar (MCS), triggered by the proton 
pulse on the spallation target at 20 Hz. The UCN counts 
were recorded for 50 ms into the time channels on the 
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FIG. 6: A MCA spectrum of the background collected with 
empty cell (red circles) is fitted to Eq.(j5} (red curve). The 
grey vertical lines indicate locations of the first 7 Bragg 
diffraction peaks from the stainless steel guide. When a mir- 
ror is used (black solid circles), the CN Bragg peaks drop in 
intensity, with some residual background scattered from the 
mirror's support frame and the vacuum chamber. UCN signal 
spectrum collected from the 4.8 K 0-D2 target is plotted (blue 
squares) and fitted to Eq.Q (blue curve). Note the relative 
intensities of Bragg peaks of the filled and the empty cell. 
The reduced Bragg peak in the filled cell is a result of the 
attenuation of cold neutrons by the presence of 0-D2 solid. 



7 



TABLE II: The expected and observed locations of the Bragg 
reflections in the background spectrum for an FCC lattice, 
a=3.6A. After re-binning, measurement wavelength resolu- 
tion is ~ O.ObA. The intensities of the [222] and [400] peaks 
were too low to be observed. 



Indices 


Expected Peak (A) 


Fitted Peak (A) 


111 


2.89 


2.90 


2 


2.50 


2.52 


2 2 


1.77 


1.78 


3 1 1 


1.51 


1.51 


2 2 2 


1.44 




4 


1.25 




4 2 


1.12 


1.14 



spectra are accumulated by the MCS for 1000 s for each 
set of experimental condition. Without the loss of proton 
pulses, 20,000 passes of spectrum are collected for each 
run. The first 1.2 ms of the detector response was dom- 
inated by the initial burst of radiation from the proton 
pulse, and is cut in the analysis. The TOF of UCN to 
travel from the production cell to the detector is much 
longer (^several seconds) than the time between cold 
neutron pulses (=50 ms), so we expect the UCN signal 
to have a uniform TOF spectrum. On the other hand, 
the measured TOF spectrum displays several prominent 
peaks, as shown in Fig. [5] These peaks originate from the 
Bragg diffraction of cold neutrons from the stainless steel 
UCN guide that intersects the cold neutron beam. The 
TOF information allows one to associate the observed 
peak to the Bragg diffraction peaks. The width of the 
peaks is due to the large angular acceptance of the UCN 
guide system. The guide was made of 316L stainless steel 
3.175 mm thick. Its lattice structure is face centered cu- 
bic with lattice parameter 3.6 A [38|. Using a total flight 
path to the detector of 24.5 m, and scattering through an 
arc of 44°, we reproduced these peaks at the correct posi- 
tions in the TOF spectrum. The expected and measured 
peak positions are summarized in Table ITT] In addition 
to these distinct peaks, there exists a long TOF tail that 
resembles the spectral shape of the incident cold neutron 
spectrum. We attribute the long TOF tail to the inco- 
herent elastic scattering of neutrons from the UCN guide. 

In order to further reduce this source of background, 
we replaced the stainless steel bend with a 45° mirror 
made of a thin nickel foil only 0.127 mm thick, housed 
inside a stainless steel tec . This reduced the amount 
of background CN by an additional factor of 4 (from 
134 mHz to 34. 8± 1.3 mHz), but it also reduced the UCN 
transmission efficiency to 73 ± 2 % of the previous con- 
figuration. Nevertheless, the signal to background ratio 
improved by more than a factor of two. The background 
TOF spectrum with the mirror is compared to that with 
the bend in Fig. [BJ 



IV. BACKGROUND SUBTRACTION 

For the data obtained in the present geometry, a two 
parameter fit to the CN-induced background and UCN 
signal is implemented. First, the model used to fit the 
raw TOF spectrum is the following: 

C(t) = C ucn + Cbg(t) 

Caen ~t~ (B CN (t)+C B ) (4) 

with 

Due to the long propagation time, the count rate of UCN, 
C ucn , should be fiat throughout the time channels. On 
the other hand, the CN background Be Nit) has detailed, 
time-dependent features due to both coherent and inco- 
herent scattering off the UCN guides. It can be modeled 
as the sum of distinct Bragg peaks (the five most promi- 
nent ones are included) and a continuous background due 
to diffuse scattering. The time-independent background 
Cb is introduced to account for the cosmic ray back- 
ground and the detector electronic noise. In the analy- 
sis, the first 10 channels are cut due to the contamination 
from the prompt radiation created by the proton pulse 
on target. Occasionally, the proton pulse was observed to 
arrive earlier than the to trigger, so we also must cut the 
final 2 channels. The positions of the diffraction Bragg 
peaks are summarized in Table HU The peaks were fit- 
ted using Lorentzian functions with different width Ti 
and coherent amplitude A| , for i = 1 — 5. The widths 
could be quite large due to the large angular acceptance 
of the UCN guide geometry. The amplitudes Af h , A mc , 
and Cb are determined simultaneously in the final stage 
of the fit using \ 2 reduction techniques. Examples of 
MCS data and fitted spectra of a typical run are shown 
in Fig. H3 

All data from background runs with an empty cell 
throughout the experiment were summed up to construct 
a parameterized background function _Bcw(i), that is 
then used as a standard background function for UCN 
production runs. Much care has been taken in subtract- 
ing this background function from the UCN production 
data. With D2 (or O2) solid in the cell, the primary 
cold neutron beam was attenuated due to elastic scatter- 
ing, reducing the population of cold neutrons which could 
elastically scatter from the UCN guide, resulting in a re- 
duced CN background. On the other hand, cosmic ray 
and electronic noise characterized by Cb presumably re- 
mains constant between runs. Therefore, the background 
function used for UCN production runs was modified to 
be: 

C bg (f,t) = ((1 - /) + fe- n ° B e c m N ptv (t) + C B , 

(6) 

where the cold-neutron associated background function, 
B^j N py (t), is derived from the empty cell data (i.e., 
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C UC n = 0), and / is the effective fill level discussed in 
sec. Ill Bl Here, / quantifies the overall attenuation of 
the CN background through the equation above, and thus 
no arbitrary scale factor is required. The intensity of the 
diffraction peaks changes between the filled and empty 
cell cases, due to the cold neutron attenuation by the 
scattering target, leading to a self-shielding effect. The 
exponential dependence characterizes the energy depen- 
dence of the cold neutron attenuation. The total cold 
neutron scattering cross-section, a, was measured from 
the in-situ cold neutron transmission measurements (as 
discussed in sec. Ill B[) . and x is the length of the cell. 

The data was fit to C(t) = C ucn + Cb g (f,t) through 
\ 2 minimization. To further increase counting statistics, 
different regions of interest (ROI) were separated and re- 
binned individually. The ROI were chosen to represent 
the distinct features of the data: the last 5 ms of the 
frame channels were mostly UCN and cosmic background 
only, whereas in the range from 25-40 ms the signal con- 
sisted of UCN and incoherent elastically scattered neu- 
trons and between 15-25 ms the signal was dominated by 
the diffracted neutrons. Six groups were fit separately 
and the final determination of the free parameters was 
made via a grid search algorithm that minimizes the com- 
bined x 2 , defined as 



E 



M eas 



CFit)iY 



(7) 



where <ii is the statistical uncertainty in the i th group. 
Typically, the reduced x 2 was close to 1, and no larger 
than 1.5. The x 2 distribution is consistent with expecta- 
tions. Nominal results of the fit are shown in Fig. [6] for 
a measurement using 0-D2 at 4.78 K. 

In the final minimization procedure, only / and C ucn 
are free parameters. A typical scan of \ 2 in the param- 
eter space is shown in Fig. which presents the anal- 
ysis of a data set collected using S-D2 at 4.78 K. The 
result shows low correlation between C ucn and /, reflect- 
ing that the UCN count was not significantly altered by 
changing the background level through adjusting the fill 
level /. The minimized x 2 is 1-26 for this temperature. 
However, if one sets the filling fraction to 1, the mini- 
mized x 2 rises to 1.70, indicating a poor fit to the data. 
The averaged filling fraction across all temperatures is 
/ = 0.61 ± 0.06 (with the uncertainty determined from 
the standard deviation across the data points at different 
temperature), consistent with the value estimated from 
the total cross-section estimate and the pressure change 
in the gas storage volume. Error bars on the values of 
C ucn and / are set by the largest extent of the 68% con- 
fidence level (A^ 2 = 1) in determining each free param- 
eter. For each 1000 s run, the UCN count rate can be 
measured to within 8% relative uncertainty. 

To facilitate the comparison of UCN production be- 
tween runs with different operational conditions, we de- 
fine the normalized UCN signal, S, as the total UCN 
count rate (C ucn determined from the fitting routine) 




3.0-10 12 



FIG. 7: A scan of x for the 4.78 K measurement. A dot 
indicates the location of the minimum, at reduced x 2 — 1-26, 
and the contours which enclose the 68% and 95% confidence 
levels area also drawn. A narrow range of UCN signal is 
returned across a wider range of filling fractions, indicating 
the relative in-sensitivity to the / parameter. 



normalized to the incident cold neutron current, Jcn, 
and the molar amount of material under study, N mo i. 



S 



1 Cucn X 512 
Nmol JcNf 



(8) 



Here Jcn is the number of total neutrons passing 
through the cell per unit time, normalized to 100 /zA of 
proton beam on the spallation target. It is related to the 
neutron flux by Jcn = A ce u<j), where A ce u is the area of 
the beam incident on the production cell and (f> is the cold 
neutron flux. Jcn is calculated from the neutron signal 
measured by the cold neutron monitor Mi, multiplied by 
the ratio of the cross-sectional area of the cell and that of 
the monitor. The filling fraction / accounts for the fact 
that not 100% of the incident cold neutrons are inter- 
secting with the target. The UCN count rate, C ucn , has 
units of count rate per channel. It is thus multiplied by 
the total channel number, 512, which corresponds to the 
whole time range of 50 ms, to calculate the total count 
rate. Normalizing the UCN count rate to the neutron 
flux monitor is a more robust method than normalizing 
to the proton beam current, as this avoids introducing 
the additional uncertainty from the inevitable varying op- 
erational conditions of the liquid hydrogen cold neutron 
moderator (such as temperature, fill level and ortho/para 
ratio). 

Since the cold neutron detectors (Mi and M 2 ) are op- 
erated in the current-integrating mode, there is an addi- 
tional step required to convert the detector signal I(t) to 
the cold neutron current Jcn passing through the cell. 
The total neutron current ranging from 1.5 to 20 meV is 
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determined by: 



>CN 



29.3 ms 

■ E 

19.5 ms 



I(ti), 



(9) 



where I{ti) is the digitized voltage signal output from the 
neutron detector in the i th time channel. The calibration 
scaling factor, 



(7.17 ±0.04) x 10 5 



neutrons 
Volt ' 



(10) 



is determined by setting the total neutron flux to the re- 
sults of a recent flux measurement in FP12 [Ic}. This up- 
dated measurement of flux using a fission chamber gives 
a total neutron flux of 2.0 x 10 7 cm _2 s _1 (with 100 [iA of 
proton current) over the energy range of interest to UCN 
production (i.e., 1.5 to 20 meV). Note that in Eq.©, 
only part of the measured TOF spectrum (between 19.5 
to 29.3 ms (i.e., 2.7 to 6.0 meV) is used to determine the 
scaling factor k. The electronic interference from adja- 
cent power lines (40 A) that supply the pulse tube refrig- 
erator compressor introduced additional noise, and thus 
the time channels with low neutron currents are excluded 
from the summation. 

A simple analysis that subtracts the background using 
the empty cell data tends to overestimate the background 
for reasons previously mentioned. However, for later runs 
in which the nickel foil mirror replaced the bend in the 
neutron guide, the background was significantly reduced, 
and the simple background subtraction routine produced 
signals in reasonable agreement with the algorithm pre- 
sented above. The simple subtraction approach uses the 
model: 



S. 



simple 



1 (£,;Ci — TY^Bi 



N, 



mol 



Jcn/ 



[11) 



where E^C, is the total number of counts measured with 
a filled cell, and the background S^S, is the total number 
of counts measured with an empty cell. As mentioned be- 
fore, the first few channels that are contaminated by the 
prompt pulse are cut. The reduction of the background 
from the self-shielding effect due to the presence of the 
UCN production medium is captured by the transmis- 
sion fraction, T, averaged over the incident cold neutron 
energy spectrum. To properly account for the fact that 
not all cold neutrons were attenuated uniformly by a cell 
partially filled to fraction /, the transmission factor is 
calculated as 



T=(l-/) + /x 



J(f>{E)dE ' 



(12) 



where T = 0.651 for the "bend" measurement (/ = 0.57) 
and T = 0.691 for the later "mirror" measurement (/ = 
0.507). 

With the proper background subtraction, we can con- 
struct a figure of merit to compare different guide ge- 
ometries. Comparison of the signal-to-noise ratio (SNR) 
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FIG. 8: UCN production in solid 0-D2 below the triple point. 
The removal of the scattered cold neutron background by the 
fitting procedure (red squares) is commensurate with a simple 
background subtraction approach (blue circles). To show the 
difference between the two data sets, the blue circles are off- 
set in temperatures by 0.25 K. The replacement of the bend 
by the mirror (white triangles) reduces the UCN transmission 
by 24 %. 



between the mirror and the bend guides can be made 
using 



SNR = 



TB' 



(13) 



For the bend SNR is 5.6 and it improves significantly to 
12.5 using the mirror. Results using the non-linear fitting 
algorithm are within the error of the simple subtraction 
method. The experimental results from both approaches 
are displayed in Fig. [8j The subtraction method is most 
useful when counting statistics for each individual TOF 
channel is so low, such that integration over all chan- 
nels is required to generate a better statistical confidence 
in the UCN signal. The agreement between the two ap- 
proaches warrants the adoption of the simple background 
subtraction method for the forthcoming analysis on S-O2 
(despite the rather unusual background shape). 

The data in Fig. [S] correspond to ortho-D2 in a 3.556 
cm long (117.4 cc) cell with bend and mirror. The solid 
was cooled over 14 hours starting from the liquid phase. 
In the liquid phase at 20 K, there was already a non- 
zero number of UCN detected well above the background 
level. Solidification begins at the triple point at 18.7 K. 
The observed UCN counts increased monotonically as the 
temperature of the production target decreased until a 
saturation was reached around 10 K. The enhanced UCN 
production was expected as the MFP of UCN increased 
with the reduction of the thermal phonon population, 
leading to a suppressed upscattering loss rate. The satu- 
ration can be understood primarily by the suppression of 
up-scattering at lower temperatures to the limit in which 
the escape length of UCN became comparable to the up- 
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The simulated UCN count rate C S i m could be pre- 
sented function of the various efficiencies through 
different parts of the apparatus: 



£elec& trans (j^source^o^J i 



(14) 



where 



is the detector efficiency, Etrans is the trans- 



port efficiency for a UCN to travel from the source 
through the guide system and finally to the detector 
active volume, n is the molecular number density, and 
Vgource is the volume of 0-D2 illuminated by the cold neu- 
trons. The cross-section a is the UCN production cross- 
section averaged over the incident cold neutron spectrum: 



1 



1 A*eV 



20 meV 



dE ucn / dE<}>(E)a{E -> E ucn ). 

Jo 

(15) 

To include the detectable VCN signal, the integration 
include neutron energy up to 1 /LteV. The transmission 
efficiency is dependent on the neutron energy. As shown 
in Fig. [TTJ some VCN with energy larger than the Fermi 
potential of the UCN guide can be present among the 
detected signals. 



FIG. 9: The temperature dependence of this experi- 
ment compared to similar experiments ( Mainz [39l]. LANSCE 
UCNA [3, PNPI PSI Hal)- All data are scaled to unity 
at 5 K. 



scattering length. In addition, cold neutron scattering 
out of the cell by the UCN source was enhanced due to 
the increased density of the cold solid. A survey of results 
from previous experiments show similar dependence over 
the same temperature range (l4[ Il6j |. with the exception 
of the data sets from LANL and Mainz (see Fig. [9J . 



V. DISCUSSION 

We built a Monte-Carlo code on the GEANT4 frame- 
work [40l | , with a complete geometry of the experiment to 
simulate the results using relevant UCN cross-sections. 
The simulation applies the physics of UCN production 
in 0-D2 , and assumes simple transport processes to track 
UCN from the source, through the UCN guide into the 
UCN detector. We simulate the expected UCN produc- 
tion using modest assumptions on the material properties 
(such as guide surface reflection specularity, UCN up- 
scattering cross-section in D2, etc.). We then compare 
the measured UCN signal S mea s to the simulated UCN 
signal S S i m , defined in Eq.© to fine-tune the guide pa- 
rameters. Furthermore, the guide parameters were con- 
strained using additional data on the guide transmission 
from the UCNA source. We assume all guide sections 
have the same surface qualities, to reduced the number 
of free parameters in the simulation. 



A. Spatial Distribution of UCN Production 

Over the energy range relevant to UCN production, the 
clastic MFP of the cold neutron beam varies from 2.8 cm 
to 5.6 cm. Since the cell length is comparable to the elas- 
tic MFP, a large probability of elastic scattering of the 
incident cold neutron beam inside the cell is expected, 
which would significantly alter the spatial distribution 
of the cold neutron flux. This distribution is modeled 
using another Monte Carlo code (MCNP5 HH). The 
code simulates the multiple scattering process of each 
cold neutron and tracks the evolution of the cold neu- 
tron flux throughout the target cell. A large fraction (~ 
50%) of cold neutrons experienced at least a single scat- 
tering, and may downscatter to UCN before escaping the 
0-D2 volume. The results of the simulation can then be 
compared to the measurements of the transmitted cold 
neutron beam collected by the detector M2 (sec III B() . 

To account for this spatial variation of the cold neutron 
flux, a position-dependent UCN production rate, R(z), is 
defined which varies as the distance, z, along the sym- 
metry axis of the cylindrical cell. To explore the position 
dependence, the target cell is sub-divided into slices and 
the production rate is averaged over the energy spectrum 
of cold neutrons in each individual slice, i.e., 

R(z) ._ 1 i: +d * dz C dE < E "> Eucn)HE, r) 

This calculation includes a divergence of the neutron 
beam of 0.5° at the guide exit (a result derived from an- 
other simulation using VITESS HI). The UCN produc- 
tion cross-section, a(E — > Ejjcn), used in the calculation 
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FIG. 10: The position dependence of the UCN production 
reaction rate, calculated from Eq. (|16|l , The integral cold neu- 
tron flux's evolution through the volume is also shown. 

is from the work of [25| , integrated over final state UCN 
energy up to 1 /j,eV. The cross-section quoted in (25| was 
scaled to the calculation using the simple Debye model 
and the incoherent approximation. The production rate 
is summed over each volume segment dV z , defined as the 
cross-sectional area of the cell multiplied by the thickness 
dz. The incident neutron spectrum, 4>{E,r), used in the 
Monte Carlo is measured by Mi. 

The position-dependent UCN production rate, R(z), 
and the cold neutron flux, <f>(z), are shown in Fig. 1101 
at Sz =1 mm resolution. Results of the simulation show 
that the cold neutron energy spectrum is not significantly 
changed by the presence of 0-D2, however, the intensity 
decreases along the length of the cell as more cold neu- 
trons are elastically scattered out of the cell. Notice that 
the spatial distribution does not follow a simple expo- 
nential decay, probably due to the fact that the cell is 
of finite volume. Finally, the cold neutron flux peaks at 
a few mm into the cell because of the locally enhanced 
cold neutron population through accumulation of elasti- 
cally scattered neutrons and the increase of the neutron 
flux in the adjacent volume slices around the scattering 
site. This spatially-dependent production rate function is 
then used to create the UCN source term in the GEANT4 
simulation. 



B. UCN Transport Efficiency 

The total transport efficiency of UCN through the ap- 
paratus can be divided into three independent compo- 
nents: efficiency of extraction from the cell, efficiency of 
transport through the guide system, and detector effi- 
ciency. The cell, cold-shield, and nickel shield lining are 
modeled with parameters provided by the known absorp- 
tion and scattering cross-sections. The up-scattering of 
UCN in these materials is assumed negligible due to their 
low temperatures and small thicknesses. The D2 volume 
is then the greatest source of loss. While the scattering 



and absorption cross-sections are well known, the up- 
scattering cross-section depends on the model employed, 
allowing for discrimination between different models of 
this effect. 

The diffusivity and additional elastic cross-sections are 
difficult to constrain. Reasonable assumptions about the 
magnitudes of the various parameters are made, and a 
comprehensive simulation is used to assess the sensitiv- 
ity of the efficiency to these parameters. To this end, a 
UCN transport Monte Carlo based on GEANT4 is used, 
which includes all of the UCN transport physics. UCN 
are created inside the cell according to the spatial dis- 
tribution above (Eq. ([To])). The initial UCN spectrum is 
proportional to v 2 dv ~ \J~EdE. Neutrons with energy up 
to 1 fieV are tracked. 

The stainless steel UCN guide walls have a neutron 
potential energy of 189 neV, and 1% diffusive surfaces. 
The loss per bounce used in our guide model is 8.5 x 1CU 5 . 
Imperfections due to welding were measured to be 2.5 cm 
long rings at the interface between the guides and conflat 
flanges, and were modeled in the Monte Carlo as 100% 
diffusive regions of the guide. The foil loss was char- 
acterized in a followup experiment to measure the guide 
transmission using the UCNA solid-D2 source. The UCN 
intensity in the guide system with and without the tee 
in place was measured, and the results are used to con- 
strain the guide parameters. The diffusivity of the foil is 
determined to be 38%, with which the simulation repro- 
duces the transport reduction of 76% when the bend is 
replaced by the nickel foil mirror. 

Finally, the loss in passing through the detector win- 
dow (a 0.51 mm Al foil) is not negligible. A transmission 
efficiency through the foil of 50% for a drop of 1 m at the 
UCNA source was measured, similar to the drop in our 
system. There exists significant additional loss beyond 
the typical absorption and elastic scattering in Al, de- 
spite the gravitational acceleration from the 1 m vertical 
drop. The diffuse clastic scattering in the Al foil is tuned 
to agree with this measurement. 

The detailed results from the simulations reveal that 
about 3/4 of the UCN produced are lost before escaping 
the cell. The loss of UCN occurs mostly on interaction 
with the cell body, which is made of Al. Aluminum was 
used to construct the cell body for its low Fermi potential 
and the ease of UCN extraction, however, its large neu- 
tron absorption cross-section is a serious source of UCN 
loss. For a UCN with velocity of 5 m/s, only 18% of 
the UCN population survives losses in the target cell and 
passage through the thin Al front window. In addition 
to the loss in the Al cell body, 14% of the UCN flux is 
lost in the solid D2 at 5 K, due mostly to upscattering. 
Upscattering loss increases to 66% at 18K. 

Furthermore, the survival of UCN through the cell 
is dependent on the UCN energy. The maximum cell 
extraction efficiency is 0.15 for UCN of 200 neV. Inte- 
grated over the UCN spectrum up to 1 /icV, the cell 
extraction efficiency at 5 K is found to be 0.07, which 
is the fraction of UCN produced which emerge from the 
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Efficiency Breakdown 




UCN energy (neV) 



FIG. 11: Top: A breakdown of efficiencies of UCN trans- 
port in the apparatus as a function of UCN energy. Bottom: 
UCN transport efficiency weighted by the production spec- 
trum f(E)dE oc y/EdE. The energy-averaged efficiency is 
obtained by integrating the curve up to the maximum energy 
of interest. The green curve is the efficiency of neutron trans- 
port through a UCN shutter used later in the experiment. 



cell, travel through the cold shield, and enter the guide. 
The transport efficiency through the guide is also energy- 
dependent. With a typical cos(f?) angular distribution 
(Eq.(2.68) in Q), the probability of transport, assum- 
ing a continuous UCN spectrum up to 1 /LteV, through 
the guide to the detector window is 0.31. Together with 
the transmission of 0.50 through the detector window, 
the 3 He(n,p) 3 H reaction probability of 0.86 and the full 
energy peak efficiency of the detector of 0.85, one esti- 
mates for an overall detection efficiency of 0.008. The 
efficiencies for UCN transport and detection are listed in 
Table. HIH 



C. UCN upscattering 

The temperature dependence of UCN production orig- 
inates solely within the production source. The UCN 
upscattering cross-section is strongly temperature depen- 
dent. The density of the target material also changes 
with temperature, but the dependence is much weaker 
by comparison. Using the upscattering cross-section cal- 
culated in our previous work(43|, our simulations did not 
reproduce the experimental data. In particular, the sim- 
ulation predicts saturation at a much lower temperature, 
around 6 K. In this calculation, the incoherent approxi- 
mation was used to treat the inelastic coherent scattering 



TABLE III: Efficiency Breakdown and UCN production Es- 
timates. Presented numbers are spectral-weighted efficiencies 
and cross-sections for UCN with with a continuous energy 
spectrum up to 1 peV to include the VCN population, and 
for UCN with energy up to the 300 neV in the guide, after 
the energy boost of 104 neV upon exiting D2. 



Parameters 


0-1 /ieV 


0-300 neV 


cell extraction 


0.07 


0.15 


guide transmission 


0.31 


0.80 


detector window transmission 


0.50 


0.50 


3 He(n,p) 3 H probability 


0.86 


0.97 


s e i ec (total energy peak) 


0.85 


0.85 


total efficiency 


0.008 


0.05 


9(b) 


1.42xl0~ 6 


1.27xl0~ 7 


No 2 o (cm _2 s _1 ) 


(44±3) x 10 30 


Expected Rate ( s - ) 


0.49±0.03 


0.28±0.02 


Measured Rate ( s _1 ) 


0.48±0.05 



in the same fashion as the incoherent scattering. 

In evaluating the cross-section for incoherent scatter- 
ing, the density of states is used to weight the contribu- 
tion of the different phonon modes. While the incoher- 
ent approximation works quite well to estimate the total 
cross-section, it might not be appropriate when evaluat- 
ing cross-sections for UCN scattering, where the phase 
space for the process is significantly limited. The lack of 
detailed Q-dependence information in the incoherent ap- 
proximation leads to significant errors in calculating the 
UCN cross-sections. 

While the incoherent approximation works quite well 
for hydrogenous neutron moderators (cr^L=80.27 b, 
CT^ C =1.7583 b), in the case of o-D 2 , the contribu- 
tion of coherent scattering is not small (er^ h =5.592 b, 
cr^ c =2.050 b). To address this concern, a new calcula- 
tion is implemented, using a full model that includes spin 
statistics and coherent inelastic scattering, as well as the 
incoherent scattering. This new calculation predicts that 
the revised UCN upscattering cross-section is a factor 
of 2-4 smaller than the calculation using the incoherent 
scattering. 

With an updated upscattering cross-section calculated 
using the full model of the dynamic structure function 
S(Q,u) for o-D 2 [45|, the Monte-Carlo simulation pro- 
duces a temperature dependence of UCN production that 
agrees better with the experimental data. Both mod- 
els, along with the background-subtracted data for the 
bend and mirror measurements, are shown in Fig. 1121 In 
Fig. [T21 to accentuate the temperature dependence, two 
sets of experimental data and the simulation results are 
scaled to agree at 4.78 K (where upscattering is small in 
both models) to eliminate guide and detector transmis- 
sion efficiency effects. At temperatures higher than 12 K, 
where the UCN upscattering cross-section is quite large 
(~ 10 b), the incoherent model is excluded at the lrj level. 
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Other effects that fead to increased clastic scattering in- 
side the source, such as the presence of cracks and voids 
leading to inhomogeneity scattering, have been studied, 
however, the higher saturation temperature can only be 
explained by reduced upscattering cross-sections. 

Increasing the elastic scattering inside the UCN pro- 
duction source (solid D2) would lead to local trapping 
of UCN and thus amplify the effects of any loss mecha- 
nism inside the source. This effect will result in a steeper 
temperature dependence on the UCN production data 
(simulated results plotted in Fig. [T2"]) . The solid D 2 was 
solidified from liquid and cooled slowly over the duration 
of 10 hours without any deliberate thermal shock. In or- 
der to quantify the elastic scattering of UCN due to the 
presence of inhomogeneity, one needs to carry out mea- 
surements of UCN MFP by systematically varying the 
cell dimension. We have done this study on the solid O2 
target, but did not perform it with solid D2. On the other 
hand, comparing the experimental data with GEANT4 
simulations with varying degrees of UCN elastic scatter- 
ing suggests that the inhomogeneity scattering of UCN 
inside the production target is no more than a few barns. 
This can be compared to the minimum elastic scattering 
cross-section, arising from incoherent elastic scattering, 
which is about 2 b. This result is consistent with the 
findings of the Los Alamos and PSI experiments on solid 
D2 sources, which are consistent with modest variations 
in the effective elastic scattering cross-section for typical 
solid D2 crystal growth. 

While the results of the simulation using the updated 
upscattering cross-section agree quite well with our data 
and the data set measured independently at PSI (shown 
in Fig.O, the data sets from the LANL and Mainz groups 
showed a much steeper temperature dependence for tem- 
peratures higher than 10 K. The difference comes from 
the different source configurations. In these two experi- 
ments, solid D2 was condensed from vapor at the end of 
the UCN guide, which was cooled below the solidification 
temperature of D2. The source was designed to reduce 
the transmission loss by eliminating a vacuum window 
which would contain the D2 . It worked quite well for low 
temperatures, however, at temperatures higher than 10 
K, the whole UCN guide was filled with D2 gas at the 
saturated vapor pressure. The additional upscattering 
from the D2 vapor could be quite large, resulting in a 
temperature dependence steeper than the simple predic- 
tion where UCN are upscattered through single phonon 
absorption. 

In this reported work, the essential focus is to use solid 
0-D2 as a source to calibrate the overall efficiency of the 
apparatus to detect UCN. In these calibration runs, the 
absolute count rate of the UCN detector can be estimated 
using Ea.([T4"|): 

Count Rate = et r ans£eiecN D2 4> a, 

where N Da = n x fV = 3 x 10 22 cm" 3 x (0.60 ± 0.03) x 
117 cm 3 = (2.2 ± 0.2) x 10 24 is the total number of D 2 
molecules in the cell. With 100 /iA proton beam current 
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FIG. 12: Results of 0-D2 UCN production experiment and 
GEANT4 simulation. Simulation parameters (guide effi- 
ciency, production, etc.) detailed in the text. 

on the spallation neutron target, the cold neutron flux 
integrated over the UCN production cross-section is 

<j) a = (2.0 ±0.1) x 10 7 cm- 2 s- 1 (1.42 x 10~ 6 b) 
= (28 ± 1.4) x lO^V 1 

Here the UCN production cross-section, a, is estimated 
following Eq.(|15p using our updated calculation, which 
includes the coherent scattering process. Note that the 
production cross-section reported in [2q is normalized 
to each atom and does not include the spin statistics 
and the rotational form factor for molecular deuterium; 
the cross-section independently reported in [44| should be 
corrected by a factor of two to properly account for the 
range of momentum transfer integration through —q ucn 
to +q UC n- With the total detection efficiency of 0.008, 
the expected neutron count rate is 0.49 ± 0.03 s _1 , which 
agrees with the measured count rate of 0.48 ± 0.05 s _1 . 
This confirms the general validity of our model of UCN 
transport. 



VI. CONCLUSION 

An apparatus to produce and measure UCN from dif- 
ferent converter materials has been constructed, the de- 
tails of which has been presented in this paper. To allow 
for useful measurements in a pulsed neutron beam line, 
a background subtraction technique has been developed 
to extract UCN production data on top of a high level 
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of background from cold neutron scattering. Using neu- 
trons on FP12 at LANSCE, the transport efficiency of 
the apparatus using 0-D2 solid as a calibration source 
has been measured. Successful application of a Monte 
Carlo model which includes detailed physics of cold neu- 
tron attenuation, UCN upscattering, and inhomogene- 
ity scattering describes the data reasonably well. With 
our detailed understanding of the apparatus, the analy- 
sis can be extended to UCN production in S-O2, where 
the production rate has never been carefully measured. 
The results of UCN production in S-O2 using this appa- 
ratus will be reported in forthcoming papers. Finally, the 
data shows an evidence for a reduced UCN upscattering 
cross-section as indicated by the higher saturation tem- 
perature. With an updated UCN upscattering calculated 
from a full coherent scattering model, the shortcomings 
of the widely-adopted incoherent approximation can be 



remedied. 
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